Emergence of quasi-metallic state in disordered 2D electron gas due to strong 

interactions 
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The interrelation between disorder and interactions in two dimensional electron liquid is studied 
beyond weak coupling perturbation theory. Strong repulsion significantly reduces the electronic 
density of states on the Fermi level. This makes the electron liquid more rigid and strongly suppresses 
elastic scattering off impurities. As a result the weak localization, although ultimately present at 
zero temperature and infinite sample size, is unobservable at experimentally accessible temperature 
at high enough densities. Therefore practically there exists a well defined metallic state. We study 
diffusion of electrons in this state and find that the diffusion pole is significantly modified due to 
"mixture" with static photons similar to the Anderson - Higgs mechanism in superconductivity. As 
a result several effects stemming from the long range nature of diffusion like the Aronov - Altshuler 
logarithmic corrections to conductivity are less pronounced. 

PACS numbers: 71.30.+h, 72.10.-d, 71.10.-w 
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I. INTRODUCTION 

The question of mutual influence of long range 
Coulomb interactions and disorder in two dimensional 
electron gas (2DEG) attracted a great attention after an 
unexpected discovery of metallic state and clear metal - 
insulator transition by Kravchenko and coworkersi^ The 
very existence of a metallic state with finite conductivity 
at zero temperature is in conflict with the weak local- 
ization theory^ which predicts that in 2D even negligi- 
ble amount of disorder localizes electrons at sufficiently 
low temperature. The theory however was firmly estab- 
lished at weak coupling or for short range interactions 
only, while the metallic state exists and the transition 
was found for rather strong coupling r s = E ee /Ep ~ 10, 
where E ee is the average interaction energy per electron 
and Ep is the Fermi energy. Therefore Coulomb interac- 
tions dominate the kinetic energy and cannot be consid- 
ered "small" . In addition to an obvious difficulty to treat 
quantitatively or even qualitatively the strong coupling, 
it is not clear which one, disorder or Coulomb interac- 
tions, should be considered as a most important cause of 
the transition to an insulating state (the corresponding 
insulating state in these cases is of " Anderson" or " Mott" 
type^). Most probably it results from a nontrivial com- 
bination of these interactions. 

The standard approach starts with a commonly ac- 
cepted argument that a long range Coulomb interaction 
after " bubble resummation" of the random phase approx- 
imation (RPA) type^ becomes effectively short range. 
Therefore one can start the treatment of disorder af- 
ter this resummation was performed. Disorder is treated 
within a similar approach in which " rainbow" diagrams 6 ^ 
"ladders and crossed ladders resummation" 7 (or, more 



systematically, the "steepest descent" approximations 8 
in the path integrals language 9 ) with interaction being 
already short ranged. In this way two kinds of mass- 
less modes determining the properties of the disordered 
electron gas are identified: diffusons (describing diffu- 
sive nature of the electron motion due to impurities) and 
Cooperons in the particle - particle channel. It is the 
last which lead to weak localization due to logarithmic 
infrared (IR) divergences in leading fluctuation contri- 
bution to conduct ivitjiiS (diffusons can also lead to IR 
divergences at yet higher orders^) 

Aa wl (T) = ^-\og[Tr/h}, 
irn 

where r is a free system relaxation time. 

More sophisticated renormalization group (RG) based 
methods using "path integral"— and "a - models"£i2ii4 
with Coulomb interactionsi^*i& were developed. Consid- 
ering high order vertex renormalization, it was found 
that there are additional logarithmic IR (Aronov - 
Altshuler 17 ) divergencies: 

Aa ee (T) = - 3/4F*) logpV/ft], 
irn 

where F* the Fermi surface average of the screened 
Coulomb interaction, leading to a conclusion that long in- 
teractions increase tendency to weak localization— This 
leads to a difficulty in understanding recent experiments 
in which apparently interactions do not necessarily lead 
to rapid increase of resistivity. Recent detailed experi- 
mental studiesi 9 ^^ clearly show that near the putative 
metal - insulator transition logarithmic terms either are 
suppressed or cancel each other (several arguments were 
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put forward in Ref. against such a fortuitous cancella- 
tion). The conductivity dependence on temperature fol- 
lows the Gold - Dolgopolov'sSi linear decrease, which at 
higher temperatures crosses over to the ballistic regime 
studied in detail recently in Ref. Generally within 
this approach the Coulomb interaction is screened first 
and the disorder effects are treated later. However recent 
electron spin resonance experiment demonstrated that 
the screening length rapidly diverges when density is re- 
duced towards the transition point. The density of states 
(DOS) at the Fermi level vanishes. It was noticed long 
time ago— that the diffusive motion of electrons slows 
down the process of screening. Therefore in the limit of 
small density and when which disorder seems to play an 
important and possibly crucial role, it is reasonable to 
start from an approximation in which the interaction is 
not rendered short range. 

With these experimental facts in mind, we reconsider 
the question of the inter-relation of disorder and Coulomb 
interactions in 2DEG within a single consistent system- 
atic approach without replacing it by a short range po- 
tential from the beginning as is done in the clean limit 
or high density. The necessarily nonperturbative ap- 
proach consists of two steps. First is a variational one 
(nonperturbative in coupling) and is similar in spirit the 
Hartree - Fock for clean metals or the BCS approxima- 
tion in superconducting metals. We find in section II 
the "best" quadratic Hamiltonian representing the sys- 
tem. On this set of Hamiltonians a quasiparticle (and 
quasi - hole) Green function is a variational parameter. 
There are possible contributions in the particle - parti- 
cle (Cooper pairs) as well as the particle - hole channels 
due to Coulomb interactions (Hartree state in direct and 
Fock state in exchange channel), while interactions with 
disorder can be treated in a similar manner with the fre- 
quency dependent relaxation time being one of the varia- 
tional parameters. Possible condensates in several chan- 
nels do not realize: of course there is no condensate in 
the Coopper channel for a repulsive interaction and there 
is also no condensation in the direct channel due to the 
charge neutrality as is shown in section IIA. However the 
strong long range exchange interaction creates (even in 
the clean case 2 ^) a dip in the DOS on the Fermi surface. 
At infinitely strong coupling the DOS on the Fermi sur- 
face approaches zero (to avoid confusion, this reduction is 
not related to the one found at higher orders for screened 
interaction in Ref. Il8l see discussion of this topic in sec- 
tion IVD). This makes the electron liquid very rigid and, 
as a result, the effects of disorder are greatly suppressed. 
This in turn leads to increase in conductivity at large 
coupling. The emergence of the above phenomenon can 
already be seen on perturbative level. The first order 
quasiparticle energy shift due to exchange is: 

S p = (p\H mt \p) =-^2 V (P - P') si S n - ^rr) ■ 

It is easily shown (section IIB) that for purely repulsive 
v(p) energy of states above the Fermi level is shifted up, 



while energy of states below the Fermi level is shifted 
down. The logarithmic vanishing of the DOS is a direct 
consequence of the long range nature of the Coulomb in- 
teraction. It is important to note that the significant re- 
duction of the DOS near the Fermi level does not mean 
that the effective mass is smaller than the band effec- 
tive mass. On the contrary, it was shown in the clean 
cas o 25 ' 26 that despite this the effective mass grows with 
coupling as was observed recently in Shubnikov - de Haas 
experiments^ We comment more on that in section IIB. 

After the variational quadratic Hamiltonian (or varia- 
tional quadratic action in the path integral formalism) is 
found, we introduce in section IIIA a systematic pertur- 
bation theory around it. In the path integral language 8,9 
it is a conventional "steepest descent" expansion with 
the variational action as a saddle point. First we intro- 
duce fields describing various possible kinds of fluctua- 
tions: diffusons, Cooperons, static photons (correspond- 
ing to the direct Coulomb interaction channel), the ex- 
change and the Cooper channel interactions. The last 
two are evidently massive as well as half of diffusons and 
Cooperons. 8 - However we find that there is a nontrivial 
mixing between photons and diffusons. The phenomenon 
is very reminiscent of the Anderson - Higgs mechanism 
in superconductivity 2 ^ in which massless Goldstone bo- 
son of phase is mixing with (dynamical) photon. As a 
result both modes become "massive". In the case of 
strongly coupled 2DEG the modes are not really mas- 
sive, the density - density correlator describing diffuson 
becoming "harder": 



uJT r + e 2 D r \p\/4Tr' 
compared to the noninteracting diffusion pole 



UJT r + D r p 2 
The (static) photon becomes RPA screened: 

1 

2|p|/e 2 + 2tt" 

The renormalized diffusion constant D r increases (non- 
perturbatively) with coupling from its noninteracting 
value of D = jir/m. Similarly the renormalized relax- 
ation time r r increases with r s . 

Equation implies that electrons at large distances 
are no longer "diffusive": they obey diffusion equation 
with first space derivative only. The Cooperon on the 
other hand still retains its typical "diffusive" pole form 
Eq. Pfl. The approximation scheme at higher orders 
therefore nontrivially combines the RPA and the disorder 
resummation on the same footing. It is important to em- 
phasize that the scheme is manifestly " gauge invariant" . 
As was shown in Refs. mil that it is very important 
to ensure gauge invariance at each stage in order not to 
miss important "vertex corrections" necessary to ensure 
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charge conservation at each order of the expansion in- 
cluding the variational stage. 

Next in section IV we turn to the "fluctuation cor- 
rections" leading to weak localization. The leading fluc- 
tuation correction to conductivity is infrared divergent. 
However, as the coupling r s grows, the correction grows 
slower than the main (Drude) contribution. This jus- 
tifies the expansion even for small values of D < 1 at 
which the standard X/D expansion is invalid (D r ^> 1 is 
however required). Therefore the crossover temperature 
at which conductivity starts approaching zero is signifi- 
cantly lower than that for non-interacting electrons. This 
temperature is estimated as a temperature at which the 
perturbation theory in fluctuations breaks down, namely 
when the correction becomes a significant fraction of the 
leading order contribution. The crossover temperature, 
according to our analysis becomes unobservably small 
since it vanishes exponentially fast with coupling (due 
to logarithmic dependence of the fluctuation correction 
on temperature serving as an IR cutoff). Therefore one 
can practically (for samples of finite albeit large size) talk 
about stable metallic state in 2DEG. Our conclusions and 
discussion of the phase diagram, as well as relations with 
other approaches are subject of the concluding section V. 



ZT Jx,t,s 



x,y,t 



r x ?TpTv(x-y)T v ?ip*p. (5) 



Here a,b = 1, ...,N r are replica indices, r is the "bare" 
relaxation time describing strength of the random po- 
tential. Let us first consider, for the sake of simplicity, 
the spin polarized case and only one "valley" (return- 
ing to the general case in section IVA), which means 
that we drop the spin indices a. The path integral for- 
mulation of the variational principle, being completely 
equivalent to the standard methods like summation of 
diagrams or Bogoliubov transformations, allows, in ad- 
dition, a convenient treatment of quantum and thermal 
fluctuations. Transforming to the Matsubara frequency 
u> n = (2n+ 1)ttT and the momentum basis for the Grass- 
mannian fields 



tit = ^ T J2J2 ex p[*(^ - 

p n 

and separating regions of phase space in which inter- 
action connects electrons near the Fermi surface, one 
obtains^: 



II. VARIATIONAL PRINCIPLE FOR 
COULOMB INTERACTIONS IN THE 
PRESENCE OF DISORDER 

A. Model and the basic approximation 

1. The model 

We consider a system of electrons with effective band 
mass m* confined to a plane interacting with each other 
and with random potential U(x): 



H 



2m' 



— fl + U{x) ) C; 



+ \ [ c\ ai c xa ^V (x - y)c\ ao c xa ^. (3) 

We set h = 1 throughout theoretical parts of the paper. 
a is a spin and valley index and v(x) is the 3D Coulomb 
interaction which has the following Fourier transform 



v(p) = --[l-S(p)] 
2e p 



(4) 



A ~ Z^Z^^pn \ "~" 1 2m* '' I '"' 



(6) 



b 



—b 



+lb ib a lb ib" 

'P—q^n t —p—r,n Yr,mYq.m 

-nT ib a ib b ] 

1 Yp—q,n Y —p—r,n Yq,m Yr,m\ ' 
fj 1 , 



+i>p- q ,ni>-p~r,n v (pH"r, m i>q, m 
+i ja p- q ^-p-r,n v (p)^ b q , m ^r,m\- 

All the fermion's momenta p, q, r are now considered to 
be around pf- Both the disorder and the Coulomb inter- 
action parts, A^ and Ac respectively, have three terms 
corresponding to direct (Hartree), exchange (Fock) elec- 
tron - hole channels and the electron - electron (Cooper) 
channel. 



Here e the dielectric constant and the last term describes 
the background ensuring charge neutrality of the system, 
v(p = 0) = 0. Passing to the standard imaginary time 
path integral formulation* and performing a well known 
replica trick^ one obtains the action: 



A[^\ = / i>. 



2m* 



2. The most general quadratic Hamiltonian and the 
Hubbard - Stratonovich fields 

A convenient way to look for the most general 
quadratic action is to perform a Hubbard-Stratonovich 
(HS) transformation introducing a field for each of the 
six channels. We should not consider the direct channel 
for disorder though since it is of higher order in number 
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of replicas N r which should approach zero (we assume 
that replica symmetry is not broken spontaneously). The 
effective action in terms of these fields is rather compli- 
cated: 



A eff = -TrLog [G^] + ^[Q^ m Q^ m + A£* ro A* 

prim 

4 E [*>~ 1 (p-p')*P'» + e^»«" 1 (p-^)eS P 'n 



6 i 
pnmJ 



pp'n 



+£=&«- 1 (p-j/)£SL n ]. 



(7) 



Here Gat is fermionic Green's function in Nambu space 
defined by 



The inverse fermionic propagator is a 2 x 2 matrix 



r -i _ I G- 1 F 



(9) 



whose elements G 1 and F are themselves matrices in 
p, n and a space. The diagonal element is: 



{p'n'\G \ pn ) — &nn'& 



ah 



P 



ft rt 

2m 

' nan _ r paO 

rp^ip—p' 'nn' u nn ,rj pp' 



+ M 1 <W + ^pn 

(10) 



where the field Q describes the diffuson, $ is a static pho- 
ton field in the direct channel, while E is an "exchange 
field" . The off - diagonal element 



(p'n' \F\pn) 



ab 



Xabp,a 

^pp'^n+n' 



(11) 



condensation of the electron - electron pairs). Therefore 
Asp = Osp = 0. This by no means indicates that there 
are no fluctuations in these channel. On the contrary 
fluctuations in the Cooperon channel play an important 
role in destroying the metallic state. 

The translation invariance in space and time (we as- 
sume that the ground state is a liquid rather than Wigncr 
crystal) and the unbroken replica symmetry (assuming 
the ground state is a disordered, possibly overcooled liq- 
uid rather then electron glass^I) implies: 



D ab 

^CprimSP 
pat 

^pp'nSP 



S S n Spp' Cp , 

5 ab 6 n S v 6. 



(12) 



We will comment on the last two (nontrivial) assump- 
tions in section V. Consequently the inverse Green's func- 
tion of fermions simplifies to 



F 



(a \n- 1 \ b \ 
\p'n' 1^* \pn I 



0, 

S ab S nn 'S( P ~p')(G^) 



(13) 



(G-r 1 = ik+-| + 



2m* 



XL 



The minimization equation for the static photon conden- 
sate d is 



(14) 



However due do neutralizing background v(p = 0) = 0, 
and the right hand side of this equation vanishes. There- 
fore <p = 0. Minimization equation for q ni e p are 



contains the A and the 9 fields describing the Cooperon 
channel in the disorder part and the Cooper channel in 
the Coulomb interaction part (if the interaction were at- 
tractive this channel would have lead to superconductiv- 
ity) respectively. This effective action should be mini- 
mized as a function of all five HS fields determining the 
fermionic Green function. The solution of this variational 
problem is discussed in the following subsection. Later 
in section HI we expand the path integral around the so- 
lution of the minimization equations to quadratic order 
(harmonic approximation) to determine elementary exci- 
tation modes and then in section IV use Feynman rules 
to compute fluctuation corrections. 



3. The saddle point equations 

Minimization of the effective action Eq. J7J poses a 
nontrivial mathematical problem. Let us first remove 
obviously irrelevant fields and functional dependencies 
using symmetry arguments. Since the Coulomb inter- 
action is purely repulsive, we assume that the electro- 
magnetic U(l ) gauge symmetry is unbroken (there is no 
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S p = -Tj2v(p- P ')G p \. 

p' ,n 



(15) 
(16) 



We start with exact solution of these equations in the 
clean case. 



B. The clean limit 

1. A major simplification in the clean limit 

In the absence of disorder r — * oo, q — * and we should 
consider the second equation (|f 6|) only. Substituting the 
Green function Eq. Qf 3p. it takes a form 



= - T E 



v(p - q) 



iuj n + /i — q 2 /2m* 



(17) 



corresponding to the diagram on Fig.^ At zero tempera- 
ture (to which we confine ourselves for most of the paper) 



p-q 




FIG. 1: The self-energy corresponding to the saddle point 
equation in the clean limit 



the summation over Matsubara frequencies results in 

Ep = -i^u(|p-q|)sign[/i- g 2 /2m* - E g ] 

i 



1 



2 (27r) a 



qv (^\/p 2 + q 2 - 2pqcos(ip) S J 



(18) 



q=0 ip=0 

sign[^ - q 2 /2m* - £,], 

where y> is an angle between fermion's momenta vectors 
p and q and p, q are their lengths. Integral over the angle 
gives 

j 



Ep — 



(27T) 



IP - ?l 



-sign[pt-g /2r 



Here if [a;] is the full elliptic integral of the first kindi^i 
The interaction with neutralizing background amounts 
to subtracting the term q = p. Now we switch to di- 
mensionless variables describing deviations of particle's 
energy from Fermi surface 



P 



v/2m*/i(l+e + ...); q= ^2m*fx(l + s' + ...), (19) 



rescaling also the variational self energy function S p = 
2fie £ . In the resulting integral equation 



K[e] 



K[-4/e 2 } 



K[e — e J sign [e 



(20) 



we extend the integration over e' to {— co,oo} and ob- 
serve that there exists a solution obeying physically rea- 
sonable property that the sign of e e is the same as that 
of e (we checked that there are no other solutions in dis- 
ordered case as well). This makes the right hand side 
independent of e e and we obtain a solution 



V2r s 



k[e'] 



V2r 



■Ki[e], 



(21) 



e'=0 



where function k\ [e] can be expressed via Meijer function, 



G 



VI 



f Ol,...,O p \ 

\ Z 6i,...,6 9 J 



31 



«i[e] 
M(e) 



r sign[e][M( £ )-M(0)]; 



(22) 



1/2,1/2,1 
0,0 



/2i. 



FIG. 2: The "rainbow" diagrams corresponding to the saddle 
point equation. The solid lines denote the fermionic propaga- 
tor, the wavy lines denote the Coulomb interactions. 



It behaves as e(log[8/e] + l)/2 at small e and as 
7rlog[2e]/2 at large e. The standard dimensionless cou- 
pling 



e 

47T6 



was introduced. 

Generally the solution of the saddle point equation cor- 
responds in terms of diagrams to summation of all the 
photonic "rainbows", see Fig. |2 However, as it is well 
known^ the sum of the diagrams except that of Fig. 
vanishes identically. This will no longer be the case in 
the disordered electron gas. The fermionic Green func- 
tion near the Fermi surface, namely at small e becomes 
nonanalytic 



G e 1 - iQ n + e(l + 



7fc l0g e } 



(23) 



and it is dominated by the interaction " corrections" . 
Here Q n — Lu n /2fj, is dimensionless frequency. 



2. Depletion around the Fermi level 

Although no energy gap (the Coulomb gap) have been 
opened within this approximation, it follows directly 
from Eq. i|23|) that the DOS right on the Fermi level 
vanishes logarithmically: 



N{e) 



Iog(l/e) 



(24) 



Therefore one can term electron gas with such proper- 
ties " very marginal Fermi liquid" . The logarithmic dip 
in the DOS is even weaker than a power normally asso- 
ciated with such a situation^ As we will see in the next 
subsection, this will naturally lead to effective reduction 
of scattering off impurities thereby increasing conductiv- 
ity and making the transition to Anderson insulator at 
least more difficult. Of course this result is "perturba- 
tive" in a sense that for the inverse propagator only one 
diagram was taken, Fig. Therefore the effect " starts" 
and can be understood at weak coupling. One can inter- 
pret the minimization equations (or the Hartree - Fock 
resummation) as a renormalization of energies of the one 
particle states due to the collective effect of the many - 
body electron - electron repulsions. As we mentioned in 
Introduction the reduction in DOS in perturbation the- 
ory can be seen from the eigenvalue shifts. Indeed the 
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TABLE I: The density of states at Fermi level JV(0) (normal- 
ized to the ideal gas DOS No = l/27r) for various couplings 
r s and diffusion constants A. 



where renormalization constant are denned by 



X 


0.1 


1 


2 


4 


8 


16 


8 


0.873 


0.395 


0.219 


0.106 


0.0476 


0.0204 


2 


0.892 


0.466 


0.268 


0.129 


0.0570 


0.0246 


0.5 


0.895 


0.557 0.344 0.167 0.0714 


0.0297 


O.f 


0.845 


0.651 


0.471 


0.241 


0.103 


0.0395 


RPA 


0.961 


0.817 


0.709 


0.623 


0.506 


0.393 



sign[e] factor in Eq. (|21|) makes it clear that energy of 
states above the Fermi level are shifted up, while energy 
of states below the Fermi level are shifted down. The log- 
arithmic singularity is a direct consequence of the long 
range nature of the Coulomb interaction. 

Obviously, if there would not be a disorder to interfere 
with the screening, the RPA type of reasoning would im- 
ply that the singularity would be "smoothed away" or 
"cured" by the quantum fluctuations corrections. Higher 
orders in coupling are increasingly singular and the sin- 
gularities should be " resumed away" A However even for 
the screened interactions there is a dip in the DOSi^ It 
cannot approach zero, but might be reduced significantly. 
We calculated N(Q) for different couplings assuming RPA 
potential instead of v(p). For several couplings the re- 
duction of DOS is given in the last line of Table I. It 
provides an indication at what degree of disorder we can 
continue to use the HF approach without encountering 
strong screening effects. As one can see, the RPA is ap- 
parently less important to the reduction of DOS already 
at r s as low as 1. 

This is consistent with results of more elaborate cal- 
culations of the renormalization constant Z defined in 
Eq. i|26[l below involving the Hubbard function in Ref.l25l 
The effect of screening is much smaller at large coupling 
due to small density of the polarizing electrons. As we 
will see in the next subsection, in the disordered case 
the " marginality" is replaced at large coupling by large 
nonperturbative renormalization of the parameters of the 
disordered Fermi liquid. 



z- v = 



z-' = 



1 + — ^Re£ re *(p,^ = 0) 
Pf op 



P=PF 



2ttN(lo = 0) ' 



(26) 



l--^ReX ret (p = p F ,cj) 



Within the Hartree - Fock approximation (without the 
RPA resummation) the retarded self energy S ret (p, oj) 
does not depend on frequency. Consequently, Zp — 1 and 
since Z < 1 for repulsive interactions 2 ^ the renormalized 
mass to* is smaller than the band mass m* (m* vanishes 
for the long range interactions, see Eq. (|24p'l. Several 
groups tried to improve this beyond Hartree - Fock using 
the RPA and the Hubbard approximation ] 25 ' 2 ^ Generally 
the wave function renormalization Zp < 1 since it repre- 
sents departure of the momentum distribution from the 
ideal Fermi gas one. Therefore the question whether the 
effective mass is larger than the band mass depends on 
which of the reduction factors spacial Z or temporary 
Zp is smaller. Ting, Lee and QuinnSi obtained finite 
monotonically increasing to* > to*, while more recent 
calculation 25 indicates that at r s < 1, to* < to* and 
become larger above r s = 1. Numerical simulations of 
the same systemSi indicate that, on the one hand side 
the renormalized mass is definitely smaller than to* at 
least for r s < 5, but, on the other hand it increases 
with r s . In recent experiments the renormalized mass 
was measured using Shubnikov - de Haas oscillations in 
magnetic fields 7 Apparently that the renormalized mass 
deduced that way monotonically increases with r s . To 
conclude, the increase of the effective mass with coupling 
does not necessarily imply that the DOS at Fermi level 
cannot drop significantly. This is important for our ap- 
proach since the drop in the DOS naively should result 
in suppression of the elastic scattering off impurities due 
to reduced phase space available. Now we return to the 
general case of disordered strongly coupled 2DEG. 



C. The saddle point equations in the general case 

1. Numerical solution 



3. Effective mass increase vs the DOS drop 

The reduction in the DOS due to the repulsive interac- 
tion is related and is sometimes confused with the issue 
of the vanishing in Hartree - Fock approximation of the 
" renormalized mass" and the Fermi liquid parameter F* 
defined by 5 ^ 



Z 

~Z~f 



1 + F[. 



(25) 



Following in the more complicated disordered case 
Eq. (|15fl - l|lt)|) the same steps as in the clean case, the 



+ -»X . £*- + 



FIG. 3: The "rainbow" diagrams involving both Coulomb in- 
teractions (wavy lines) and interactions with disorder (dashed 
lines) in the general case 
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FIG. 4: The dependence of e e /e on e for various values of r s 
and A. The solid lines are the numerical results, the dashed 
lines are the results of r a expansion, and the dots are the 
results of 1/A expansion. 



FIG. 5: The dependence of g„ on uo for various values of r 3 
and A. The solid lines are the numerical results, the dashed 
lines are the results of r a expansion, and the dots are the 
results of 1/A expansion. 



equations for the scaled (dimensionless) quantities (E p = 
2fie £ , q n = 2ny/rq UJ , to = 2/icD) at zero temperature are 



1 



(2ir)2\J £ ( £ + e £ ) 2 + (£ + ^) 2 ' 

e' + e E 



V2tt 2 



(27) 
r(28) 



where A = /it and function k[e] was defined in Eq. (|20H . 
The approximation corresponds to summing up the 
whole set of "rainbow" diagrams involving both Coulomb 
interactions and interactions with disorder, Fig. El 

The results for the real part of self energy for A = 
0.5, r s — 1 (where in the non-interacting case Anderson 
localization is very effective) and A = 8 , r s = 0.01, 1, and 
16 (still a quasi-metal with small localization effects) are 
given on Fig.0J The corresponding imaginary part of self 
energy is given on Fig. |SJ We solved the equations for 
various values of coupling r s between weak coupling up to 
r s = 24 (experimentally the metal - insulator transition 
is observed from r s — 15, 1 to r s — 25 in GaAs/AlGaAsS 
heterojunctions and around r s = 10 in Si MOSFETS*^). 
Solid lines are simulation results. The dashed lines are 
results of expansion in r s briefly described in Appendix, 
while dots are the results of the next to leading order in 
1/A also summarized in Appendix. 



At zero frequency the large coupling value of q^^o is 
much smaller than the non-interacting result 1/(4A). We 
define the renormalized relaxation time via 



2qu. 



4Ag' 



(29) 



One observes on Fig. [5] that for r s < 10 the disorder 
parameter q u is almost independent of lj in the whole 
region where it is important for calculations of integrals 
over Green's functions (namely when q < uj). Therefore 
one can still consider the system as a disordered liquid 
without temporary dispersion. It is a good approxima- 
tion to simplify the analysis by considering a variational 
principle with constant q u = q. 

The saddle point equations in that case take a simpler 
form (after integration over u>) 



1 = 



1 



1 



(2tt)2A J e ( £ + e e f+q?' 
V2r 



(30) 



7T q 
k\s — £'|sign|e'| | — — arctan— 

2 \£ "T Cr 



(31) 

From Fig. 0| we observe that at very large coupling it is 
similar to the "very marginal Fermi liquid" of the clean 
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TABLE II: Renormalization of the inverse relaxation time q 
for various couplings r s and bare diffusion constants A. 

















X 


0.1 


1 


2 


4 


8 


16 


8 


0.028 


0.013 


0.0070 


0.0034 


0.0015 


0.00067 


4 


0.057 


0.028 


0.015 


0.0075 


0.0033 


0.0015 


2 


0.12 


0.056 


0.034 


0.017 


0.0073 


0.0032 


0.5 


0.47 


0.29 


0.18 


0.086 


0.037 


0.015 


0.1 


2.4 


1.6 


1.3 


0.67 


0.27 


0.10 



0.08 




-0.08 -0.04 0.04 0.08 



co 



case. The interaction "correction" E £ dominates and is 
not proportional to e, therefore it does not reduces to 
a mere renormalization of the density of states on the 
Fermi level. At finite, but large coupling the renormal- 
ization of the density of states is still large, see Table I. 
We calculated it as explained in section IID, Eq. (J2HJ). 
At large coupling the quantity diverges very fast. The 
limiting value of q at zero frequency for various " bare" 
diffusion constants and couplings are given in Table II. 

Before trying to further exploit the solutions of the 
saddle point equations, we would like to explicitly show 
that even perturbatively the reduction in the DOS due 
to exchange can be clearly seen. 



D. Reduction in density of states 



0.5 
0.4 


1 i 1 

" - -v . A,=8 
r =l 


1 

(b) 


0.2 






0.1 











-0.08 -0.04 

CO 


0.04 0.08 



FIG. 6: The DOS N(u>) (normalized to the ideal gas DOS 
No = l/27r) for (a) fixed r s = 8 and various A, (b) fixed A = 8 
and various coupling r s . 



The fact that at small e (deviation of fermion's energy 
from n) and large coupling r s , S £ 3> 2fis leads to signif- 
icant rearrangement of the DOS in the vicinity of Fermi 
surface. The DOS N(u) is related to the retarded Green 
function by: 



N(ui) 



Im / G: e >). 



(32) 



For frequency independent q u it reduces to 



N{lo) 



n 

1 

2^2 



Il m y 

7T * J !i) — 



1 



lo — 2/iEp — S p + iq 



Q 2 



(33) 



The DOS at lu = for various couplings are given on 
Table I, while the DOS as a function of u> for fixed r s = 8 
and various A on Fig. 6a and for fixed A = 8 and various 
coupling r s on Fig. 6b. At large coupling the renormaliza- 
tion of limiting value of r r is very large and q approaches 
zero. The integrand in Eq. (|33|l becomes delta function 
and the DOS is determined by the derivative of e e with 
respect to e at e = 0: 



N(w) 



1 1 



2tt 1 



u)— £+e e 



(34) 



Using the simplified saddle point equation Eq. (fST]) the 

derivative of e F is: 



de 



e=0 



k'[e] 



e>0 



arctan • 



(35) 

The asymptotic of the function n'[e] is as follows. At 
small e it is negative and divergent, k'[e] w — l/2s. As 
we have seen in the clean case this lead to vanishing DOS. 
However in the presence of disorder the divergence is cut 
off since the second multiplier vanishes at small e as e s /q. 
This in particular means that in the case of disorder the 
RPA improvement is not necessary as long as the cutoff 
due to disorder is larger than the cutoff due to screening 
(see discussion in subsection IIB). 

At large e the integral converges rapidly since k'[e] « 
— 7r/2e 2 . This can be seen perturbatively as well using 
results of the previous subsection. To leading order in r s 



2\/2r| 
7T 2 qi°) 



k' [e]e 



(i) 



£>0 



where cf ^ and ejp are given in Appendix Eqs. (| A. 1|> - 
121. The frequency dependent correction to the DOS 
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(compared to ideal gas Nq = l/27r) is perturbatively 



6N(u) = N(lo) - N Q = 



(2D - e)ei 1} 



8tt 2 A 2 

In particular on the Fermi surface 
SN (0) = 



(2D - ef + qi°P 



1 

where P, p are a total and relative momenta of the two 
electron state respectively. The relation to indices of the 
© fields used in effective action Eq. Q is obvious 

P=Pi+P2] p = pi-p 2 - 



where cf- 1 ^ is given in Appendix Eq. I|A.3|) . The results 
for the DOS are given in Table I. The effect of reduction 
of the DOS at the Fermi level due to repulsion increases 
fast with increasing r s and depends weaker on A. As we 
mentioned, in the clean case the DOS approaches zero 
at any coupling no matter how small, which is clearly an 
artifact of neglecting screening at this stage. However the 
disorder effectively "averages" the distribution of states 
making it finite. At the limit of large disorder the DOS 
should approach the ideal gas one. Repulsive interaction 
works to reduce the DOS near the Fermi surface and is 
expected to make scattering by impurities less effective. 
We will see in the next section that e' is related to the 
renormalized value of the diffusion constant: D r = A/e'. 



III. A SYSTEMATIC EXPANSION IN 
FLUCTUATIONS AROUND THE SADDLE 
POINT. 

A. Classification of fluctuations around the ground 

state. 

All the Hubbard - Stratonovich fields correspond to 
elementary (harmonic) excitations of the system. The 
spectrum of these excitation is determined by a quadratic 
term in expansion of the effective action Eq. JJJ around 
the saddle point determined by a solution of the min- 
imization equation Eqs. H27() - (|28() . This quadratic form 
should be diagonalized to find the spectrum. We first sort 
out evidently massive modes which do not play a role in 
subsequent discussion. These are exchange E and the 
Cooper exchange O fields. In the next two subsections 
we address photons, diffusons and Cooperons. 

Since the ground state does not break spontaneously 
the electric charge U(l) symmetry the charged fields O 
and Cooperon A do not mix with neutral fields E, dif- 
fuson Q and static photon <£. In this sector the inverse 
propagator elements (second functional derivatives of the 
effective action with respect to relevant fields) is: 



pi p . 



D 



Hi-) 



Pp 



n m 



D 



AH 



Pp 



= 5{P-P')5 nn ' S(p - p'Hp)- 1 

fi(P P )Gp/2 + p/2 Gp/2-p/2 

+S( P + p')G 
= 5{P - p')s n '+ m '' n G n ' 



i 



The 09 element of the inverse fluctuations propaga- 
tor indicated that in this channel there are no massless 
modes. This is evident at small coupling since the diag- 
onal first term dominates, but since the Coulomb inter- 
action is repulsive is true for any coupling. The mixing 
with Cooperon cannot turn it to a massless mode. We 
therefore discard the field O in what follows. Now we 
turn to the neutral fields. 

Second derivative of the effective action with respect 
to E is 



P'p' 



D 



EE 



,,,,, = 8(P-P')5 nn '(6(p-i/)v(p)- 1 
- 5(p-p')G P/2+p/2 G P/2 _ p/2 ).(37) 



It is again a massive mode and its mixing with other 
neutral field is unimportant. More interesting are the 
diffuson Q and the static photon <E> fields: 



n m 
p> 



D 



QQ 



a. in 
V 



jn— m,n' — m 



'S(p-p') 



/m /~vrr 



D. 



<I'Q 



nm 
P 



D 



= -6 



5 u '5{p-p') 



l.n—m 



S(p-p') 



Erin s-im 

1 



(38) 



(39) 



v{p)~ 



T 



p-\-q q 



(40) 

To conclude there are three fields which might have po- 
tentially massless modes. While Cooperon cannot mix 
with the other two, diffusons and photons can. Even if 
both the diffuson and the photon are massless, after mix- 
ing the massless modes could turn massive. We study this 
phenomenon next. 



B. Anderson - Higgs mechanism for diffusons 

1. Matrix elements of the photon - diffuson inverse 
propagator matrix 



At small frequency and momenta around Fermi surface 
p/2+p/2 1 - r p/2-p/2' we will use the asymptotic expressions for the solution of 
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the saddle point equations 



qui ~ qu=o = Z^Vrq = 

S e w = 2^(Z~ 1 - l)e, 



(41) 



where renormalizations of the "disorder efficiency" 1/t 
and of the inverse density of states at Fermi level Z^ 1 
were introduced in Eq. (|26|) . As can be seen from Table I 
and Table II at large coupling they are quite large. 

Then one computes standard diagrams in Eq. (|38() for 
matrix elements of the inverse propagators involving dif- 
fusons and static photons at certain fixed momentum p: 



Ann 



= S 



n — m,n —m 



- - 2_, ^p+q^q 



cn~m,7i — m rcnn 



[S r ' 

—6{—nm)B{p 1 n — m)], 



(42) 



where at small to and p, it contains well studied func- 
tion "bubble" integral having the following asymptotic 
at small frequencies and momenta: 



B(p,l) 



q 



LU t T r 



D rP 2 . (43) 



for n < 0, n + 1 > 0. Within an approximation of "renor- 
malization" Eq. (|41|) one obtains: 

D r = 

Beyond this approximation the value of effective diffu- 
sion constant can be calculated from the values given in 
Tables I and II. As usual, we discard the massive same 

frequency sign ++ and modest and concentrate on 

different frequency sign excitations. 

It is convenient to this end to rescale the photon field 



%n = VV(P)T%U- (44) 

The mixing and the photon inverse propagators are: 



2. Eigenvalues and eigenmodes: physical photon and 
diffuson 

The Q$ inverse propagator matrix is "blocked" for dif- 
ferent photon frequencies I = n — m (we take n > 
m < 0), namely we can consider a single value of I. For 
a fixed frequency I > the range of possible n is limited 
to — I < n < and the (Z + 1) x (Z + 1) matrix has the 
following form: 



/ a 
a 

\b b 



b 



(48) 



where 

a = 1-B(p,l),b 



^/vfrjfB(p,l),c=l-v(p)TL(p,l). 

( I'M 

Eigenvalues of this matrix are: the (Z — 1) times degener- 
ate a ~ UJiT r + D r p 2 (original "massless" diffusons before 
mixing) and two nondegenerate eigenvalues 



A± = - [a + c ± y/(a - c) 2 + 4ZZ> 2 



corresponding to eigenvectors {1, 1, 1...1, a±} with a± — 
(A± — a) jb. Their asymptotic at small momenta and fre- 
quency is: 



A+ « l + «(p)/(27r), 



A_ « 2nv{p) 1 (uiT r + D r p 2 ) + D r p 2 



(50) 



The form of the A+ eigenvalue means that physical pho- 
ton mode <&p = Yl n =-l Qp + a +^p propagator (rescaling 
back by yJv{p)T) is the RPA photon propagator exhibit- 
ing Debye screening: 



D. 



<I>Q 



D. 



,i„i. 



= -S l ' n - m S(p-p')^v(p)Te(-nm)B{p, n-m 

(45) 

=S u '6(p-p')[l -v(p)L(p,l)}, (46) 



), 



— u>— -u) 

^p ^-p 



v~ x {p) + 27r' 



(51) 



The second is a "symmetric" in n superposition of dif- 
fuson modes = J2n=-i Qp + a_$p is no longer the 
usual diffusion pole: 



where the Lindhard function^ has an asymptotic behav- 
iors: 

L(p, I) = T^B{p, n, n - Z) « (-1 + ujir r + ...)/(2ir). 

n 

(47) 

One notices in Eq. 1)45(1 that, due to the factor 0(—nm), it 
is precisely the massless different frequency sign diffusons 
that mix with photon. The mixing strength is large since 
it is determined by the same bubble integral that appears 
in the diffuson's inverse propagator Eq. (|42|) . Now we will 
invert this matrix and find its eigenmodes. 



Qp Q—p 



1 



ujT r + D r p 2 + D r v(p)p 2 /(2n) 

1 



ujT r + D r e 2 p/ (47r) 



(52) 
(53) 



It becomes harder (less singular) than in the standard 
treatment in which the mixing between photon and dif- 
fuson is neglected. It is interesting to note that in the 
standard treatment at small couplings the mixing is not 
neglected as far as photon's propagator is concerned. One 
uses the RPA propagator Eq. I]5ip. despite the fact that 
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(a) (b) (c) 

FIG. 7: The Feynman rules for (a) fermion-fermion-diffuson, 
(b) fermion-fermion-Cooperon, (c) fermion-fermion-photon 
vertex. The solid and wavy lines denote the fermion and pho- 
ton propagator respectively, the double solid lines with oppo- 
site (same) direction arrows denote the diffuson (Cooperon) 
propagator. 



Propagators for relevant modes supplemented by ver- 
tices, the fermion - fermion - diffuson 

r- --- 

and the fermion - fermion - photon <E> 

i/^* = — * 

constitute Feynman rules shown on Fig. 7(a) and (c) 
respectively. However fields $ and Q do not correspond 
to " modes" or " bosonic excitations" of the system due to 
mixing between them discussed in detail in the previous 
subsection. One can still use these fields in calculation 
considering them as a vector and their propagator as a 
matrix. 



the same off diagonal matrix element Eq. H39JI is simulta- 
neously responsible for the essential modification of the 
diffuson. 

The inverse of a matrix of the type of Eq. (|48|l generally 

is: 



1 



lb 2 



-b 



V 



-b -b ... a J 

Matrix elements of the propagators therefore are: 

{n,n + l\D QQ \n',n' + 1) =P 1 5 nn , + P 2 ^-, (54) 

ZTTl 



(n,n + l\D Q $\l) = P 3 



(55) 



C. The Cooperon propagator 

As we mentioned in section IIIA due to its charge the 
Cooperon does not mix with photon or any other neu- 
tral field. It is massless for different sign frequencies and 
massive for the same sign frequencies. The strong cou- 
pling however influences its propagator beyond the ev- 
ident renormalization A — » D r , r — > r r . Substituting 
the expression of the fermion propagators Eq. (|13fl into 
Eq. I|36|) the inverse propagator of the excitation is: 

(rom|D^^« \nm) = 1 — 9(—nm)B(p, n — m). (61) 

Numerical solution of the saddle point equations substi- 
tuted into Eq. (|36[) show that the dependence is quadratic 
only till certain momentum at which it saturates, see 
Fig. El for r s — 1,2 and 4. This is of importance later 
when we estimate the quantum correction to conductiv- 
ity in section IV. 



(l\D^\l) = P 4 , 



(56) 



where functions P and their asymptotic at small fre- 
quency and momentum are: 



Pi 
Pi 
P.; 
Pi 



1 _ 1 

a T r oj + D r p 2 
lb 2 



v(p) 



a(ac — lb 2 ) T r uj + D r p 2 + D r p 2 v(p) 

\Thdb \J lov (p) 

ac — lb 2 T r uj + D r p 2 + D r p 2 v{p) ' 

a T r uj + D r p 2 



lb 2 r r uj + D r p 2 + D r p 2 v{p) ' 



(57) 
(58) 
(59) 
(60) 



We will see in subsection HID that the diagonal part of 
the diffuson describes density fluctuations and therefore 
mixing with photon makes electrons nondiffusive at large 
densities. 




FIG. 8: The dependence of B(p,u = 0) on p for various r s 
and A. The solid lines are the original values of B(p,ui — 
0), while the dashed lines are their corresponding quadratic 
fitting values. 
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TABLE III: Derivative B' /tau 2 of the Lindhard function in the presence of disorder for various values of r s and A 





0.1 


1 


2 


4 


8 


16 


8 


12.2 


3.13 x 10 2 


3.33 x 10 3 


6.06 x 10 4 


1.40 x 10 6 


3.12 x 10 7 


4 


5.78 


1.14 x 10 2 


1.15 x 10 3 


2.11 x 10 4 


5.14 x 10 5 


1.29 x 10 7 


2 


2.74 


40.4 


3.77 x 10 2 


7.02 x 10 3 


1.81 x 10 5 


4.93 x 10 6 


0.5 


0.620 


4.75 


34.8 


6.35 x 10 2 


1.90 x 10 4 


6.18 x 10 5 


0.1 


0.117 


0.605 


1.79 


25.3 


9.15 x 10 2 


4.10 x 10 4 



Note that the excitation remains massless even at 
strong coupling. This follows from very general consid- 
erations. Consider the "bubble" diagram 

B (^i,p) = -7-2 / G qn G q+v<m , (62) 

(ZTT) Jq 

where I — m ~ n. At zero momentum using the saddle 
point equation Eq. 115(1 one obtains: 

If 1 

B{u hP =0) = / (G qn - G qm )— ■ 

(27T) Jq I N + 9m - q n \ 



(63) 



This approaches 1 in the limit of zero frequency as long 
as there is a jump in from negative to positive fre- 



quencies. The propagator of the different frequency sign 
Cooperons at small frequencies and momenta is: 



D AA , 



1 



2quj + B'p 2 ' 

where B' denotes a derivative of the bubble integral 
dB(0,p) 



(64) 



B' = 



dp 2 



p=0 



(65) 



Assuming quadratic momentum dependence of B(p,uj), 
direct calculation leads to 



B' = 



1 



2^(2^)2 J [$2 + ( £ + e ^)2] 



/ \2 



(66) 



J 



The values of B' for various couplings r s and bare dif- 
fusion constants A determining the Drude conductivity 
are given in Table III. The perturbative expression for 
this quantity is 



B' 
where 
C = 



-Xt 1 



1+rJ - 12\<fV 



2 5/2 A 3 



c 



(67) 



(4A) 
+2G 



((4A) 1 



3/2,3/2 
1,1 



r>23 ( 1 



r 43 



V(4A) 2 



2,2,2 

3/2,3/2,3 



(4A) 2 



2,2,2 
3/2,3/2,4 / ' 



(68) 



( ;: : Za ) oftheMeijer 



functions^ It rises fast with r s , while being weakly 
dependent on A. Finally the fermion - fermion - Cooperon 
vertex is 



l/>7/jA 



and it completes the Feynman rules shown on Fig. 7(b). 



D. Density - density correlator and conductivity to 
leading order 

1. Modification of diffusive motion due to strong 
interaction 



One of the most important characteristics of 2DEG 
is the density - density correlator describing the diffu- 
sive nature of the charge carrier's motion in a disordered 
medium. It is closely related to dielectric function and 
polarizability of 2DEG. The correlator is given in Mat- 
subara formalism by 



11 ,92 



p+q 



First we use the Feynman rules stated above to calcu- 
late the density - density correlator at the leading or- 
der. In the limit of small frequencies the contributions 
come from diagrams (a)-(d) on Fig. EI They are —L(p, I), 

-u t B(p, 2 Pi(l + P4)/2n, 2y/u; l v(p)/^L(P, l)B(j>, l)P 3 , 
and —v(p)L(p,l) 2 P4, respectively and can be combined 
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to 




(b) 

GO 





(c) 



(d) 



FIG. 9: Contributions to the density-density correlator at the 
leading order. 



into 



Here expressions for L, B and P4 are given in Eqs. (j47jl . 
(I62jl . (j57jl - lj60jl and the " noninteracting" correlator is de- 
fined by 

( r , vi B(p,l) 2 

Its asymptotic at small u> and p is: 

D r p 2 



X(^;,p) 



LOlTr + D rP 2 + D r p 2 v{p) / '(2tt) ' 



Therefore not surprisingly it is proportional to prop- 
agator of the "diagonal" diffuson defined in Eq. 1)52(1. 
The diffusive behavior dominates short range fluctuations 
only on scale smaller than s = 2e/e 2 . On a larger scale 
the last term in the denominator is linear in p and is 
therefore larger than the standard diffusion term. This 
makes diffusion less long range although in 2D it does 
not become a short range one. The scale was introduced 
by Si and Varma^fi and we will comment on connection 
to their work in section V. 



2. The leading (Drude) contribution conductivity 

The DC conductivity can be read off the density - den- 
sity correlator using the relation 



lim lim — - 



-X{u,p) 



u>->o 2tt (1-B(w,0)) 2 

Z7T 



(69) 



which follows from the Kubo formula^ Here we used 
the asymptotic of the bubble integral B Eq. (|47|l and the 



imaginary part of self energy q and derivative of the bub- 
ble integral B' defined in Eqs. (|15fl and l|65j) respectively. 
Results are given in Table IV for various r s and A. One 
observes that at large coupling the " Drude" conductivity 
increases considerably compared with the noninteracting 
one. We explain this by reduction of interaction with dis- 
order (g is much smaller than its noninteracting value of 
1/2t) despite the reduction in density of states (B' larger 
than its noninteracting value of — fi). At small r s using 
Eqs. HA.2j l- lA.3jl of Appendix and Eq. trT7j) one obtains 



p 2 / 95/2 \3 

CT= e —\(l-±r s \q (1) +r s ^^C 



where C is given in Eq. I|68jl . The leading order contribu- 
tion dominates at small couplings and disorder. However, 
the theory has zero modes - Cooperons. Therefore possi- 
ble IR divergencies might render the leading order results 
invalid at large coupling or disorder. In principle for zero 
temperature and infinite samples the results are invalid 
for all couplings. Our next task is find a range of param- 
eters and temperature (or sample sizes) in which the IR 
divergencies at the next order are still small compared to 
the main contribution. 



IV. SUPPRESSION OF WEAK LOCALIZATION 
BY THE LONG RANGE INTERACTION 
EFFECTS 

A. The saddle point expansion and the spin-singlet 
approximation 

In this section we describe some of the corrections 
around the variational ground state found in section II 
and used in section III to calculate several physical quan- 
tities. The steepest descent expansion in terms of Feyn- 
man diagrams is quite standardA^ We briefly describe it 
introducing N s identical " spin" components to show that 
the expansion might be interpreted as "1/N S " expansion 
(spin might include other degeneracies like multiple val- 
leys in Si). The action including the spin indices a is 
given in Eq. JSJ. It is a peculiar feature of the disor- 
dered Coulomb problem that the leading order in 1/N S 
vanishes for two entirely unrelated reasons. The direct 



TABLE IV: The Drude conductivity a/2ive 2 for various values 
of r a and A. 



X 


0.1 1 


2 


4 


8 


16 


8 


9.85 51.4 


1.68 x 10 2 


7.15 x 10 2 


3.39 x 10 3 


1.53 x 10 4 


4 


4.79 22.0 


70.3 


3.01 x 10 2 


1.47 x 10 3 


7.18 x 10 3 


2 


2.33 9.31 


28.7 


1.24 x 10 2 


6.22 x 10 2 


3.20 x 10 3 


0.5 


0.553 1.61 


4.44 


19.0 


1.03 x 10 2 


5.77 x 10 2 


0.1 


0.107 0.246 


0.470 


1.79 


10.5 


68.3 
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contribution in the disorder part vanishes due to the fact 
that it is of higher (second) order in replicas N r as well: 



/ j '<TlCT 1 '- J "<T 2 <T2 



The direct A 2 contribution to the Coulomb part 

/ EE G ^ x W x - (*> y) = 

Jx >V a o 1 o 2 

vanishes due to neutralizing background J v(x — y) =0 
(and under assumption of homogeneity). Therefore lead- 
ing terms are of order N s . The free theory action is also 
of order N s . Therefore all the terms in action are of 
the order N s and it plays a role of the "loop expansion 
parameter" and comes always in combination with 1/h. 
This Hartree - Fock logic is not rigorous however. There 
is an assumption involved: it is assumed that all the Hub- 
bard - Stratonovich fields which in general are tensorial 
are dominated by their singlet part: 



Q" 



2 Q. 



In this paper we will make such an assumption. Therefore 
we neglect, for example, triplet channels in the physical 
case of N s = 2. 

The partition function (suppressing for simplicity the 
replica indices and writing explicitly just one of the HS 
fields) or any observable is expanded around the saddle 
point 



Z = I e NsA "" [Q] 



« f eiq>[N s A eff {Qsp} + ^QD- 1 [Qsp}Q + AA[Q]}, 

where A^4[Q] contains all the cubic, quartic and higher 
order terms in Q. From this Feynman rules are read 
and they scale compared to Fig. in the following way: 
HS fields' propagators are proportional to N s „ fermion 
loop also has N s , while fermion - fermion - boson vertex 
is 1/y/Ng The leading order contribution of conductivity 
considered so far is of order A s and we will consider order 
A? = 1 in the next section. 



B. Fluctuation correction to the density - density 
correlator and conductivity 

1. Density - density correlator 

The correction to density - density correlator at two 
loop order which contributes to the small frequency limit 
(the only ones needed for subsequent calculation of the 
conductivity) is given on Fig. ITUI 




FIG. 10: The correction to the density- density correlator at 
two loop order. 



(n, n + l\D QQ \m, m + l) G™£ G™+ r G^+ l +r 
(n, m'\D^^ r \m + l,m' + I) 
(m', m + 1\Dqq\u', n +1) . 

All the other diagrams are regular as u> — > , hence 
they do not give contributions to the DC conductivity. 
Near the Fermi surface one "disentangles" the momenta 
flowing in the central loop, see Fig. 
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(l-B(p,l)) 



■(p.oEr 



B(r,l) : 



where 



B*(p,l) 



E/~m s-in-\-l r^n r^n-\-l 



(70) 



(71) 



The integral over r is logarithmically infrared divergent in 
2D and, as usual jiS signals breakdown of naive perturba- 
tion theory and appearance of weak localization effects. 
We assume an IR cutoff (to be defined more explicitly 
below) and will use this expression to calculate conduc- 
tivity. 



2. Weak localization 

The fluctuation correction to conductivity using Kubo 
formula is: 



8a = lim 



B(0, If 



c^o 2tt (1-5(0,0)" 



r^E— I 



B(r,iy 



where the derivative of B4 Eq. I|71ll is defined by 



B' 4 



dB 4 {p,l = 0) 



d{p 2 ) 



p=0 



After some algebra it takes a form 
1 f (e + e £ ) 2 



B 4 = 



M^) A Je (? + (£ + e e yy 



-(i + 4) 2 



^E E ■ 

q,r nmn'm' 



Values of the coefficient B' 4 for various couplings and dis- 



-n(n + l))9(—n'(n' + l))B(p, 2 order strength are given in Table V. 
I 
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TABLE V: The values of B' 4 /t 4 for various r a and A 



V- 


0.1 


1 


2 


4 


8 


16 


8 


3.77 4.64 x 10 2 


1.60 x 10 4 


1.25 x 10 6 


1.47 x 10 s 


2.11 x 10 10 


4 


3.46 


2.86 x 10 2 


9.05 x 10 3 


7.21 x 10 5 


8.92 x 10 7 


1.28 x 10 10 


2 


3.18 


1.70 x 10 2 


4.78 x 10 3 


3.85 x 10 5 


5.14 x 10 7 


7.78 x 10 9 


0.5 


2.73 


53.6 


1.04 x 10 3 


8.10 x 10 4 


1.35 x 10 7 


2.58 x 10 9 


0.1 


2.44 


27.4 


1.23 x 10 2 


6.65 x 10 3 


1.51 x 10 6 


4.72 x 10 8 



Before discussing physical implications of the correc- 
tion we provide perturbative result for B' A . 



B i = 2Ar 4 



2 5/2 A 2 



l + r s \(-2 2 5^ + t— 
\ 6n 



where 

C 4 = 



° r<22 
2?A 4 33 



((4A) ; 



3/2,3/2 
1.1 



/ 1 



12G : 



(4A) 5 



2,2,2 

3/2,3/2,4 



Q/^23 



(4A) 2 



2,2,2 
3/2,3/2,3 

2,2,2 
3/2,3/2,5 



The quantity is increasing very fast with coupling r s and 
decreases slowly with disorder strength. 
Returning to conductivity one obtains 



5a 



e 2 2g 2 B' i 

27T 7TT 



rdr 



1-B(r,l) 



In 2D the integral is dominated by small momenta. 
Therefore at very low temperature we can use approx- 
imation of Eq. 



6a = 



e 2 2q 2 B\ f puv rdr _ e 2 q 2 B' 4 p 2 uv 



2?r tttB' 



_ = _ 1 " 4 log i 

r 2 2ir tvB> & p 2 R 



(72) 



As usualifi it is cut off in both infrared and ultraviolet. 
The infrared cutoff for the weak localization logarithmic 
divergence can be set by finite temperature 



P 2 ir 



2irTm*/T 



2( q /^)(B'/r 2 ) 
or finite size L of the sample 



Pir 



The ultraviolet cutoff i: 



vlv 



,30 



2m* 



2tt 

T' 



i{h,1/t}. 



(73) 



C. Crossover temperature 

Let us find a temperature at which the perturbation 
theory in "loops" or 1/A^ S breaks down. At this temper- 
ature the Drudc conductivity is significantly reduced by 
fluctuations and one conservatively estimates it as set- 
tling of the weak localization (the Anderson insulator) 
regime. It is estimated by equating leading (Eq. (JEU) 
and the fluctuation correction (Eq. (|72|I % 1 to conductiv- 
ity times a factor R of order 1 at finite temperature (or 
sample size): 



a Q = —4N s q 2 B'/T = R6a 

Z7T 



e 2 q 2 B' 
2tt tttB 



7 log 



2m* 



I 



2nT w im* / t 
2(q/^)(B'/r 2 ) 



(7,4) 



where we used the large disorder value in Eq. (|73|l . There- 
fore 



T 



2qB' 

7rr 3/2 



with a dominant argument of the exponential being 

_ AttN s (B') 2 
RB' 4 

Considering first the bare diffusion constant as fixed one 
observes that as r s increases the temperature first rises 
due to preexponcntial factor, but then after reaching a 
maximum exponentially drops at large r s . On the other 
hand fixing r s the temperature quickly drops as A in- 
creases. In experiment what is usually varied is den- 
sity of electrons. In this case as density gets lower both 
the diffusion constant becomes smaller and r s becomes 
larger. The overall effect is that for clean samples and rel- 
atively large r s the quasimetallic state is stable to very 
long temperatures due to reduction of the DOS at the 
Fermi level. Note that the trajectory in the A, r s space of 
the experimental setup is itself dependent on the DOS. 40 
This complicates the actual comparison since effects of 
screening cannot be neglected in experiments to date, as 
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FIG. 11: Vertex correction due to interactions with disorder 
(dashed lines) in the conventional notations 18 (first line) and 
in the present paper notations (second line). 

was discussed in section II. Qualitatively however the pic- 
ture is that at large coupling the metallic state survives 
effect of scattering off impurities due to the reduction in 
the DOS. 



D. Higher order effects. Aronov - Altshuler effect 
revisited. 

1. Higher order corrections to the vertex function and 
conductivity 

It was shown in Ref. that in perturbation theory 
that when one sums up all the corrections to the vertex 
part (the ijnpP condensate where p is the density field 
that couples to static photon) shown on Fig. El it be- 
comes proportional to the diffusive pole 1/(wt + Dp 2 ). 
The same expressions are also shown in Fig. 1111 in our 
notations as a sum of leading and the next to leading 
terms in the steepest descent expansion. In this pic- 
ture however the diffuson propagator is considered in the 
noninteracting theory. The vertex part enters high-order 
diagrams creating logarithmically divergent corrections 
which strengthen (in the singlet sector) weak localiza- 
tion. The major diagrams involving the singular vertex 
part contributing to conductivity are given on Fig. ^] 
(some other contributions cancel, see Ref. [HI). The phys- 
ical interpretation of this phenomenon is that electrons 
scatter coherently on Friedel oscillations due to density 
fluctuations^ 

Without the crucial pole factor the Aronov - Alt- 
shuler's corrections do not diverge in the infrared. In 
the strong coupling regime considered in this paper, due 
to mixing of diffusons with static photons, the vertex 
is given in Fig. 1 181 has much softened small momentum 
asymptotic: l/(wr -I- Dp 2 ). It can be easily seen that 
this softening is quite enough to render the contributions 




FIG. 12: Major contributions to the conductivity in pertur- 
bation theory 

like those on Fig. 1121 (which is of the order 1/N S namely 
higher than the weak localization one) finite. One there- 
fore would ask how this can be understood diagrammat- 
ically in terms of conventional disorder coupling^ ? The 
point is that the mixing effectively sums up diagrams 
to all orders in Coulomb coupling r s Fig. ^5] Each one 
of these is divergent, while their sum is not. This is 
quite analogous to disappearance of IR divergencies due 
to long range photon "chains" after the RPA diagrams 
are summed. 



2. Density of states near Fermi level at large coupling 

The Aronov - Altshuler's corrections to conductivity 
are directly related to the downturn cusp in DOS due to 
Coulomb interaction. In 2D the cusp is given byl& 

5N(e) oc log[e] 

again due to the renormalization of the vertex, see 
Fig. E] It was claimed that this is precisely what was 
observed in tunnelling junction experiments in disor- 
dered metal films^i^ In our approach, due to Ander- 
son - Higgs mechanism, this renormalization is greatly 
reduced. An alternative explanation at very strong cou- 
pling and significant disorder might be the leading order 
reduction of DOS discussed in section IID, see Fig. 
This does not contradicts the experiments in metals since 




FIG. 13: Vertex corrections due to mixing of diffusons with 
static photons in the present paper notations (first line) and 
in the conventional notations 6 (second line) 
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in these experiments disorder is large (even very large), 
while the coupling r s is quite small. Of course if the 
density is sufficiently high the screening can no longer 
be neglected and the Aronov - Altshuler's effect becomes 
dominant. However, as we mentioned before, in very 
clean 2DEG samples the disorder can reduce the screen- 
ing and our approach of neglecting the screening at the 
leading order becomes more appropriate. In this case 
higher orders will not be large enough to undermine this 
assumption. One will get again a reduction in the DOS, 
but for entirely different reason. A related issue is emer- 
gence of the Coulomb gap commented on in the next 
section. 



V. SUMMARY AND DISCUSSION 

To summarize we present a consistent gauge invariant 
approach to disordered strongly interacting electron gas 
in 2D. Physically the basic phenomenon is the reduction 
in the DOS at Fermi level due to strong Coulomb repul- 
sion. This in turn suppresses both screening and scat- 
tering of impurities stabilizing the metallic state against 
weak localization effects. Formally the approach consists 
of two steps. The first is variational (or "self consis- 
tent"): the most general quadratic states were consid- 
ered and one with minimal energy identified. At this 
stage the "RPA" screening is neglected assuming it is 
sufficiently weakened by disorder so that higher orders 
are small. The second step is steepest descent pertur- 
bative expansion (which also can be identified as an ex- 
pansion in parameter 1 /N s with N s number of spin com- 
ponents or valleys). Although the general philosophy of 
the steepest descent expansion is not drastically different 
from the one adopted in other works (see for example 
Ref. |3), two rather independent observations were made. 
The first is that the exchange part of Coulomb interac- 
tion leads at strong coupling to a significant reduction 
of the DOS near the Fermi surface and to via the re- 
duction suppresses the disorder effects. The second is 
that when the steepest descent expansion procedure, if 
followed consistently, mixing between static photons and 
diffusons not only causes Debye screening of the pho- 
ton, but leads in addition to a softening of the diffusion 
pole. This in turn leads to a number of observable conse- 
quences like a significant modification of the vertex part 
and consequently of the Aronov - Altshuler contribution 
to conductivity. The contribution becomes regular in the 
strong coupling, not logarithmic (the absence of the effect 
of the DOS reduction due to interaction in Ref. 17 can 
be traced to an approximate calculation of the exchange 
diagram in their Eq. (16) ). 

In this section we discuss several general questions and 
assumptions and relation of our work to other attempts 
to incorporate the long range Coulomb interactions into 
the theory of disordered electron gas. 

The theory of disordered electron gas relies to large ex- 
tent on existence of massless collective modes, diffusons 



and Cooperons. It is tempting to interpret these exci- 
tations as Goldstone bosons of some symmetry breaking 
(with several complications arising from the quenched 
disorder, see Ref. E^l . The a model approach initi- 
ated by Wegner— and othersiii4 long ago and devel- 
oped and applied to the Coulomb interaction case re- 
cently by Baranov et al^£ start from and assumption 
that the G = Sp(2N) <£> Sp(2N) symmetry of free dis- 
ordered electron gas is spontaneously broken down to 
diagonal subgroup H = Sp(2N), where N enumerates 
replica, spin and Matsubara indices. We have shown in 
section III however that diffusons mix with photon and 
become "harder" than standard Goldstone bosons. This 
might signal that the a - model approach should be modi- 
fied to incorporate the Anderson - Higgs mechanism. Ac- 
tually the need for such a modification can be found in 
recent remarkable work of Ref. and we comment on 
this now. 

Unfortunately the presence of strong Coulomb inter- 
actions explicitly breaks a subgroup of G. An exam- 
ple of the explicitly broken symmetry transformation is 
Sippn = ip- p n i 1 ^p.n = — V'-p.-n f° r positive Matsubara 
frequencies n > and Sip pn = i/>_ Pi _„, Stp p n = — -0_ P: _ n 
for negative Matsubara frequencies n < O.The symmetry 
is broken by both the frequency term ^ ip pn (—iuj n ) i\) pn 

and by the Coulomb interaction term. However while the 
breaking by the frequency term is "soft" and insignif- 
icant, as far as static quantities like DC conductivity 
are concerned, it was showni& that the Coulomb inter- 
action effectively represented on the a - model level by 
the "square of trace" operator (Eq. (2.1) in Ref. Ilfll) is 
relevant and cannot be reduced to a soft breaking. Bara- 
nov et al. notice that at large distances the diffusion is 
suppressed which coincides with our Eq. (O- At short 
distance scales the electrons are diffusive. We believe 
that the Anderson - Higgs mechanism can be treated ap- 
proximately within the a - model approach as long as the 
mixing is small. This is well known problem in Quantum 
Field Theory^ under the name of "gauged a - models". 
These issues might become clearer when the present ap- 
proach is extended beyond 2D (say to 2+e). Work on this 
is in progress. A related issue is understanding the dif- 
ference between diffusons and Cooperons. Within the er- 
model approach the Sp(2N) (g) Sp(2N) symmetry forces 
the Cooperon and the diffuson propagators to be same. 
It is precisely an explicit (not spontaneous) symmetry 
breaking due to Coulomb interactions that makes the 
hardening of the diffuson (mixing with photon), while 
leaving Cooperon intact possible. 

The fact that Coulomb interaction modifies diffusion 
at large distances was also discussed in Ref. |s3 and this 
might lead also to suppression of weak localization. The 
work however was criticizedji£i22i^ that the vertex parts 
were not taken into account or alternatively the treat- 
ment is not gauge invariant. Our work explicitly shows 
that despite the fact that diffuson is "harder" (although 
still massless) at large distances, the Cooperon is not. 
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Therefore although weak localization is suppressed, the 
suppression is much weaker and completely different. 
The logarithmically divergent contribution to conductiv- 
ity includes Cooperon. 

It was shown by Efros and Shklovskii^ on basis of 
a heuristic argument with plausible assumptions about 
the nature of the localized electronic states (neglecting 
their overlaps) that there should be a Coulomb gap in 
the strongly interacting electron gas. As in the case of a 
models with Coulomb interactions^ it is not clear from 
the first orders in our scheme whether the reduction of the 
density of state is along the line of their argument. We 
believe this is unlikely due to the fact that they neglect 
the effect of exchange on the "states" ipi defined there. 
It is also not clear at this point whether the opening of 
Coulomb gap that Baranov et deduce on the basis of 
the 2 + e expansion is related to the reduction in density 
of states due to exchange. 

One can extend the approach presented here to the 
"self consistent" scheme initiated by Vollhardt and de- 
veloped to include Coulomb interactions by Sadovsky. 39 
This will allow quantitative study of the insulating state 
and of the Coulomb gap. Note however that the "gap" 
equations of Ref. |3^ follow the perturbative Aronov - 
Altshuler contributions, while the "self consistent" form 
of our conductivity contributions Eq. (|FT5|l and Eq. ff^ 
will contain different diagrams. The work on this is in 
progress. 

At last wc briefly comment on two general assump- 
tions made. The first is the spatial homogeneity. It is 
clear— that clean and even disordered 2DEG— at suf- 
ficiently strong coupling become inhomogeneous Wigner 
crystals or "glass". It was even speculated 2 ^ that the 
Wigner crystallization (which occurs around r s = 40 in 
clean systems) might be related to the observed metal - 
insulator transition. Following general argument can be 
advanced against such a scenario. It has been observed 
recently that in several clean systems of thermally fluc- 
tuating repelling objects the homogeneous state (liquid 
or gas) exists down to zero temperature. One such a sys- 
tem is one component classical plasma— Another is a 
system of vortex lines in type II superconductors— The 
latter is quite analogous to 2DEG. The difference is that 
thermal fluctuations should be replaced by quantum and 
bosonic field by fermionic (statistics is quite unimportant 
in the low density limit though). To be sure the energy 
of the solid is lower, so below the melting point the liq- 
uid state is metastable (in conventional liquids for which 
in addition to repulsive interaction there is a long range 
attractive force, the metastable state ceases to exist at 
spinodal point). It is reasonable to assume (and it was 
demonstrated recently^) that disorder favors homoge- 
neous state over a structured crystal). Therefore tran- 
sition to a Wigner crystal or glass state would occur at 
much higher couplings than metal - insulator transition 
and the relevant state is homogeneous as was assumed in 
the present paper. 

Another assumption commonly made is that the 



replica symmetry used to derive our starting point Eq. (0 
was assumed to be unbroken. This means that we ne- 
glected a possibility of " electron glass" — This is a distant 
possibility in the quasimetallic state since, as we argued 
in the paper, reduction in the DOS due to long range in- 
teractions make disorder less favored. Eventually in the 
insulating state glassy behaviors will eventually prevail. 
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APPENDIX: APPROXIMATE SOLUTION OF 
THE MINIMIZATION EQUATIONS 

1. Expansion in small r 3 

From Eq. (|31|) we observe that at small coupling 
r s <C 1, e e starts from the first order: e s = 

sign[e] (r s e^ + Substituting this into Eq. I|3U|) im- 
mediately gives the leading term for q u : 



g2» =signM- = signMg< ). (A.l) 



Therefore we expand q^ — sign[w] + r s qu' + 

Then the leading order contribution to e £ can be com- 
puted: 



(i) - yl 



K\e — e 



e' + e e , 



v Ju>o,e> - {e' + e e ,) z + {Q + qi°)y 
= J KF-£jsign(e ) I - - Arc tan — 

(A.2) 

We use here the constant q^ approximation. Substituting 
this into Eq. I|3U|) one obtains: 

^ = _3V2(g< >) 2 f k[e] 



2tt 2 J e (e 2 + 4g( )2) 
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■ 215/2^3 A 2 33 



(4A) 2 



3/2,3/2 



(A.3) 



Z P9 



/ ai,...,a p \ 



are defined 



where the Meijer function G£ 

in Ref. |^. We observe that it is negative, namely the 
long range interaction reduced the effect of disorder! Di- 
agrammatically the minimization equations sum all the 
rainbows including both disorder and interaction, Fig. 2. 
In perturbation theory one evaluated a diagram with one 
photonic rainbow and arbitrary number of disorder lines, 
namely diagram Fig.l with disordered Green's function. 
Its imaginary part is precisely The actual expan- 

sion parameter is ^/2r s /n 2 rather than r s as can be seen 
from comparison of the perturbative and exact solutions. 
Therefore perturbation theory breaks down completely 
at r s ~ 10. 



expand in 1/A. The leading order is the clean solution 



already discussed in subsection IIB: e. 



[o] _ jo] 



4 1] /a- 



with e s = \/2r s Ki [s]/tt. The leading term for q u 



.is: 



4tt 



(A.4) 



The correction to dispersion relation can also be com- 
puted in a quite compact form: 



n[e — e'] — k[e + e'} 



Jo] 



Jo] 



(A.5) 



2. Expansion in small 1/A 

It is important also to obtain an analytical solutions 
of the minimization equations at large r s . This turns out 
to be possible for relatively clean case in which we can 



The actual expansion parameter is l/47rA rather than A 
as can be seen from comparison of the perturbative and 
exact solutions. Therefore perturbation theory breaks 
at A ~ 0.1. The results are marked by dotted lines on 
Fig.0|and|SJ Generally numerical results agree with both 
perturbative expansions. 
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